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Abstract 



We derive upper bounds for the smallest zero and lower bounds for the largest zero of Laguerre, Jacobi and 
Gegenbauer polynomials. Our approach uses mixed three term recurrence relations satisfied by polynomials 
' corresponding to different parameter(s) within the same classical family. We prove that interlacing properties of 
, the zeros impose restrictions on the possible location of common zeros of the polynomials involved and deduce 
strict bounds for the extreme zeros of polynomials belonging to each of these three classical families. We show 
numerically that the bounds generated by our method improve known lower (upper) bounds for the largest 
(smallest) zeros of polynomials in these families, notably in the case of Jacobi and Gegenbauer polynomials. 
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1. Introduction 

The behavior of the zeros of orthogonal polynomials has attracted significant interest from both theoreticians 
and numerical analysts since the first results were proved by Markov and Stieltjes in the late 19th century. A 
•'-H . wide range of tools and techniques have been developed to analyse different properties of the zeros including 
' Markov's theorem on the monotonicity of zeros (cf. [17]); Sturm's comparison theorem for the zeros of solutions 
^ . of second order differential equations (cf. [22]); Obrechkoff's theorem on Descartes' rule of sign (cf. [20]); and 
the Wall- Wetzel theorem on eigenvalues of Jacobi matrices (cf. [25]). Discussion of results and techniques 
related to the properties of the zeros of both general orthogonal polynomials and classical families of orthogonal 
polynomials can be found in [6] and [18]. 

It is well known that if {pn}J^o orthogonal sequence, then the zeros of p„ are real and simple and each 

open interval with endpoints at successive zeros of p„ contains exactly one zero of Pn-i] a property called 
the interlacing of zeros. Stieltjes (cf. [22], Theorem 3.3.3) extended this interlacing property by proving that 
if m < n — 1, provided Pm and p„ have no common zeros, there exist m open intervals, with endpoints at 
successive zeros of p„, each of which contains exactly one zero of Beardon (cf. [2], Theorem 5) proved 
that one can say more, namely, for each m < n — 1, if Pm and p„ are co-prime, there exists a real polynomial 
Sn-m-i of degree n — m — 1 whose real simple zeros together with those of pm, interlace with the zeros of p„. 
The polynomials Sn~m~i are the dual polynomials introduced by de Boor and Saff in [5] or equivalently, the 
associated polynomials analysed by Vinet and Zhedanov in [24]. Related results are proved by Segura in [21]. 

In recent years, authors including Ismail and Muldoon (cf. [14]), Krasikov (cf. [16]), Gupta and Muldoon (cf. 
[12]), Ismail and Li (cf. [15]), and Dimitrov and Nikolov (cf. [7]) have developed interesting methods including 
the use of chain sequences and the derivation of inequalities for real-root polynomials to refine and improve 
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upper and lower bounds for extreme zeros of classical orthogonal polynomials. In this paper, we prove that if 
two polynomials and Pn+i, k < n, satisfy a three term recurrence relation of a certain type, common zeros 
of Qk and Pn+i cannot occur at consecutive zeros, nor at the largest or smallest zero, of Pn+i- We show how 
this knowledge of the allowable location of common zeros, together with the Sticltjcs intc^rlacing property, leads 
to strict lower (upper) bounds for the largest (smallest) zero of p^+i in any orthogonal sequence {pn}^o- We 
apply our results to three families of classical orthogonal sequences, namely, to sequences of Jacobi, Gegenbauer 
and Laguerre polynomials and compare our bounds to those derived by other authors for these classical families. 

2. Results 

The first result that we use is proved in [10]. 

Theorem 2.1. (cf. [10, Lemma 3.2]) Let {pn}^o ^''^V sequence of polynomials orthogonal on the (finite or 
infinite) interval (c, d). Let Qn-i be any polynomial of degree n — 1 that satisfies 



for for each fixed n € N where is a constant (depending on n) and a„ and bn are functions with 6„(x) 
for X e (c, d) . Then, for each fixed n e N, 

(i) the zeros of gn-i are all real and simple and, together with the point An, they interlace with the zeros of 
Pn+i if gn-i andpn+i are co-prime; 

(a) if g-n-i and Pn+i are not co-prime, then 

a) they have one common zero located at x = A„; 

h) the n—1 zeros of gn-i interlace with the n remaining (non-common) zeros ofpn+i', 
c) if Wn+i < ■ ■ ■ < wi denote the zeros of Pn+i, then An = Wk for some k G {2, . . . , n}. 

Corollary 2.2. For each fixed n ^'H, An is a strict lower bound for the largest zero of Pn+i and a strict upper 

bound for the smallest zero of pn+i- 

ReniEirk. Equation (1) and the conditions of Theorem 2.1 hold when gn-i = Pn-i by virtue of the three term 
recurrence relation satisfied by any orthogonal sequence {pn}^^o and hence Corollary 2.2 is a generalisation of 
Theorem 2 in [3] while Thciorcm 2.1 (ii) is a generalisation of Theorem 3 in [2] since it accommodates the case 

when Pn-i and Pn+i have a common zero. 

Our main result is a generalisation of Theorem 2.1. 

Theorem 2.3. Let {pn}^=o be any sequence of polynomials orthogonal on the (finite or infinite) interval {c,d). 
Let gn-k be any polynomial of degree n — k that satisfies, for each n G N, and k S {1, . . . , n — 1}, 



where f{x) ^ for x G (c, d) and H{x), Gk{x) are polynomials with deg(Gk) = k. Then, for each fixed n gN, 

(i) the n real, simple zeros of Gkgn-k interlace with the zeros of Pn+i if gn-k andpn+i are co-prime; 

(ii) if gn-k and pn+i are not co-prime and have r common zeros, counting multiplicity, then 

a) r < min {k, n — k} ; 

b) these r common zeros are simple zeros of Gk; 

c) no two successive zeros of pn+i , nor its largest or smallest zero, can also be zeros of gn-k', 

d) the n — 2r zeros of Gkgn-k, none of which is also a zero ofpn+i, together with the r common zeros of 
gn-k and Pn+i, interlace with the n + 1 — r remaining (non-common) zeros ofpn+i- 



9n-l{x) = an{x)pn+l{x) - {x - An)bn{x)pn{x) 



(1) 



f{x)gn-k{x) = Gk{x)pn{x) + H{x)pn+l{x) 



(2) 



Proof of Theorem 2.3. Let Wn+i < ■ ■ ■ <wi denote the zeros of Pn+i- 



(i) Prom (2), provided Pn+i{x) ^ 0, we have 



Pn+l{x) 



H{x) + 



Pn+l{x) 



(3) 



Pn+l{x) j^^X-Wj 

where Aj > for each j e {1, . . . , n + 1} (cf. [22, Theorem 3.3.5]) 




Therefore (3) can be written as 



f(x)gn-k{x) 

Pn+l{x) 



n+1 



Gk{x)Aj 



H{x) + Y^ 



X — Wj 



Xj^Wj. 



(4) 



Since Pn+i and Pn are always co-prime while Pn+i and gn-k are co-prime by assumption, it follows from 
(2) that Gk{wj) 7^ for any j G {1,2, . . . ,n + 1}. Suppose that Gk does not change sign in the interval 
Ij — {wj-^-i , Wj) where i G {1, 2, .... n}. Since Aj > and the polynomial H is bounded on Ij while the 
right hand side of (4) takes arbitrarily large positive and negative values on Ij, it follows that gn-k must 
have an odd number of zeros in each interval in which Gk does not change sign. Since Gk is of degree k, 
there arc at least n — k intervals [tVj^i, Wj), j G {1, . . . ,n} in which Gk does not change sign and so each 
of these intervals must contain exactly one of the n — k real, simple zeros of gn-k- We deduce that the 
k zeros of Gk are real and simple and, together with the zeros of gn-k, interlace with the n + 1 zeros of 

Pn+l- 

(ii) If r is the total number of common zeros of and gn-k counting multiplicity then each of these r zeros 
is a simple zero of Pn+i and it follows from (2) that any common zeros of gn-k and Pn+i must also be zeros 
of Gk since p„ and Pn+i are co-prime. Therefore, r < min{fc, n — k} and there must be at least (n — 2r) 
open intervals of the form Ij ~ {wj+i,Wj) with endpoints at successive zeros of Pn+i where neither Wj+i 
nor Wj is a zero of gn-k or Gk{x). 

If Gk does not change sign in an interval Ij = {wj+i,Wj), it follows from (4), since Aj > and H is 
bounded while the right hand side takes arbitrarily large positive and negative values for x G Ij, that 
gn-k must have an odd number of zeros in that interval. Therefore, in at least (n — 2r) intervals Ij either 
gn-k or Gk, but not both, must have an odd number of zeros counting multiplicity. On the other hand, 
gn-k and Gk have at most {n — k — r) and (fc — r) real zeros respectively that are not zeros of Pn+i- We 
deduce that there must be at most (n — 2r) intervals Ij = (wj+i, Wj) with endpoints at successive zeros 
Wj+i and Wj of Pn+i neither of which is a zero of gn-k- It is straightforward to check that if the number 
of intervals Ij = {wj+-i,Wj) with endpoints at successive zeros of Pn+i neither of which is a zero of gn-k 
is equal to n — 2r, this is only possible if no two consecutive zeros of Pn+i, nor the largest or smallest zero 
of Pn+i, can be common zero(s) oi Pn+i and gn-k- This proves a) to c) and d) follows from c). 



Corollary 2.4. For each fixed n G N and each fc G {1, . . . , n — 1}, the largest zero of Gk is a strict lower bound 
for the largest zero ofpn+i- Similarly, the smallest zero of Gk is a strict upper bound for the smallest zero of 

Pn+l- 



1. Equation (2) and the conditions of Theorem 2.3 hold when gn-k = Pn-k- This can be seen (cf. [2, Thm 
4]) by iterating the three term recurrence relation satisfied by any orthogonal sequence {pn}'^=o- Theorem 
2.3 (ii) is more general than Theorem 5 in [2] for orthogonal polynomials since it caters for the case when 
Pn-k and Pn+l have common zeros. 



Remarks. 



2. Gibson (cf. [11]) proved that if {pn}^=o is any orthogonal sequence and m < n — 1, then p„ and pm can 
have at most niin{m, n — m — 1} common zeros and he showed that this upper bound is sharp. Theorem 
2.3 (ii) a) generahses Gibson's result. 

3. It turns out that common zeros of two polynomials is a frequent occurrence within most classical orthogonal 
families (cf. [1]) and also occurs when the two polynomials arc members of different orthogonal sequences 
within the same family (cf. [8], [9] and [10]). It is therefore important to accommodate the possibility of 
common zeros in the analysis of Stieltjes interlacing of zeros of polynomials. 

3. Applications to classical families and comparison with existing bounds 

Markov's Theorem proves that the zeros of classical orthogonal polynomials vary monotonically with an increase 
(or decrease) in the parameter(s). Mixed three term recurrence relations of the type given in (1) and (2) provide 
functional equations linking polynomials from the same classical orthogonal family but corresponding to different 
values of the paramcter(s). From the perspective of deriving good boimds for the extreme zeros, those mixed 
three term recurrence relations that involve the largest possible parameter difference between the polynomials, 
or alternatively, no parameter difference but the largest possible difference in the degree, while still preserving 
interlacing properties for their zeros, are particularly useful. 

In the discussion below, we shall denote the zeros of the polynomial Pn+i of degree n + 1 by Wn+i < ■ ■ ■ < wi. 

(i) (a) The sequence of Laguerre polynomials {i"}^o orthogonal on the interval (0, oo) with respect to 
the weight function a;"e~^ for a > —1. Upper bounds for the smallest zero Wn+i of the Laguerre 
polynomial , a > — 1 , namely 

(a + l)(a + 2) 

< n + a + 2 

and 

(a + l)(a + 3) 

< 2n + a + 3 

were obtained by Hahn (cf. [13]) and Szcgo (cf. [22, cqn. (6.31.12)]) respectively. Gupta and 
Muldoon recovered these bounds in [12, eqns. (2.9) and (2.10)] and the bounds also follow from 
Corollary 2.2 and Theorem 1 (iv) and (v) in [9]. 

Replacing a by a + 1, n by n — 1 in [9, eqn. (13)] we have 

x^L'^tlix) = ia + 'S){n + a + l){{a + 2){a + 4)- {2n + a + 2)x)L'^tl{x) 

+n((n - l)x^ + 2{a + 3)(n - l)x - (a + 2){a + 3)(a + 4))L'^+\x) 

and using [9, eqn. (4)] together with the three term recurrence relation for Laguerre polynomials (cf. 
[22]) we obtain 

x^L'^tlix) = {n + a + l){{a + l)4-{a + 2)2{3n + 2a + 5)x + {n + a + 2)2x'^)L^{x) (7) 
-(n + 1) ((a + 1)4 - (a + 2)2(3n + a + l)x + {n- l)nx'^) L^+i{x) 

where 

(a)fe = a(a + 1) . . . (a + fc - 1), fc > 1, /c e N 
is Pochhammer's symbol. The smallest zero of the quadratic polynomial coefficient of LJJ in (7) is 

(a + 2)2(3n + 2a + 5) - ^y{a + 2)2(9(0 + 2)2 + 2(2a + 5)(a2 + 5a + 10)n + (5^2 + 25a + 38)n2) 

2(n + a + 2)2 

and this is a sharper upper bound for the smallest zero of L^^^ than the upper bound 



(8) 



(g + l)(a + 2)(a + 4)(2n + a + 3) 
^ (a + l)2(a + 2) + (5a + ll)(n + l)(n + a + 2) ^ ^ 

derived by Gupta and Muldoon in [12, eqn. (2.11)] as illustrated in table 1. 



Table 1: Comparison of upper bounds for the smallest zero of L"_^^(x) for different values of n and a. 



Values of a and n 


a = 340.56, n = 12 


a = 65.3, n = 20 


a = -0.9, n = 100 


a = —0.9, n = 3 


Smallest zero of L'^^i 


251.815653 


27.67770854 


0.00103831 


0.0259078 


New bound (8) 


275.856 


32.6035 


0.00103832 


0.0259079 


Bound (9) given in [12] 


309.594 


36.6163 


0.00103833 


0.0259082 


Bound (6) given in [12] 


319.258 


41.8125 


0.00103909 


0.0259259 


Bound (5) given in [12] 


330 


51.111 


0.00108803 


0.268293 



Remark. As observed in [12] , it is evident from (5), (6) and (9) that the upper bounds for the 
smallest zero will be sharpest for a close to —1. 

(b) The lower bound for the largest zero wi of L'^_^-^{x), a > — 1 

wi > 3n - 1 (10) 

is given by Neumann in [19] while 

wi > 4n + 4 + a - 16^271 + 2 (11) 

in [4] is due to Bottema and 

wi > 2n + a + l (12) 

was obtained by Szego (cf. [22, eqn. (6.2.14)]). Note that (12) also follows from Theorem 1 (i) in [9] 
and Corollary 2.2. 

By iterating the three term recurrence relation for Laguerre polynomials we obtain 

(a + n- l)(a + n)L^_2(x) 
= (a;^-2(2n + Q!)a; + 3n^ + 3an + a^ + l)L^(a;)-(n+l)(2n + a-l-a;)L^+i(a;) (13) 

and the largest zero of the quadratic polynomial coefficient of in (13) yields the improved lower 
bound 

wi > 2n + Q + \/-n? +an + l (14) 

for the largest zero of 

Remark. It is interesting to note that although (14) is an improvement on previous bounds for 
most values of n and a, the bound (10) due to Neumann, which does not depend on a, compares 
favourably when both n and a are small. 



Table 2: Comparison of lower bounds for the largest zero of different values of n and a. 



Values of a and n 


a = 340.56, n = 12 


a = 65.3, n = 20 


a = -0.9, n = 100 


a = -0.9, n = 4 


Largest zero of L^^^ 


469.74362252 


172.7701885 


591.362 


11.1263 


New bound (14) 


429.612 


146.616 


298.654 


10.7606 


Bound (12) given in [22] 


365.56 


106.3 


200.1 


8.1 


Bound (11) given in [4] 


310.976 


45.6081 


175.69 


-31.4964 


Bound (10) given in [19] 


35 


59 


299 


11 



(ii) The sequence of Jacobi polynomials {Pn'^}i^=Q is orthogonal on the interval (—1, 1) with respect to the 
weight function (1 — a;)"(l + .t)^ for a, /3 > —1. 

It was proved in [10, Thm 2.1(c)] that for any fixed n G N and a, /3 > —1, (1) holds with 



,a+4,/3 , _ 2n(n + Q + /3 + 3) + (a + 3)(/3-Q) _ ^ 

' 2n(n + a + /3 + 3) + (a + 3)(a + /3 + 2) Pn+i - ^n+i- 



It follows from Corollary 2.2 that 



2(a+ 1)(q + 3) 

2n{n + a + P + 5) + {a + :i){a + P + 2) ^ ' 

Similarly, from [10, Cor 2.2(c)], (1) holds with 

9n-i - ^n-i , - 2n(n + a + ^ + 3) + (^ + 3)(a + ^ + 2) ^'''^ ^"+' " 
and we deduce that 



^ _T , 2(/3 + l)(^ + 3) 

'^"+^ ^ ^+2n(n + a + /3 + 3) + (/3 + 3)(a + ^ + 2) ^""^ 

Remark. It should be noted that for any fixed n £ N, (15) is a particularly sharp lower bound for 
the largest zero of P^+i when a is close to —1 and similarly (16) provides a sharp upper bound for the 
smallest zero of P^^i when /3 is close to —1. Asymptotically in n, (15) and (16) yield wi > 1 — 0(;^) and 
Wn+i < -1 + 0(;jj-) respectively. 

In table 3 we compare (15) to the bounds given by Szego in [22, eqns. (6.2.11) and (6.2.12)] 

2n + /3-a 

> 2n + a + P + 2 ^''^ 

wi > — — for (3 > a. (18) 

n + a + 1 



Table 3: Comparison of lower bounds for the largest zero of P"^i (x) for different values of n, a and 13. 



Value of n 


n = 3 


n = 15 


n = 11 


n = 18 


n = 3 


n = 3 


Value of a 


a = -0.9 


a = -0.9 


a = 30.9 


a = 30.9 


a = 30.9 


a = 0.9 


Value of P 


/3 = -0.8 


P = -0.8 


P = -0.8 


6 = 32.8 


P = 32.8 


P = 132.8 


Largest zero of P^^i 


0.984119 


0.999143 


0.141417 


0.71025 


0.30295 


0.989891 


New bound (15) 


0.984109 


0.999142 


-0.0507338 


0.590098 


0.182432 


0.989162 


Bound (17) given in [22] 


0.968254 


0.993399 


-0.179298 


0.372665 


0.110181 


0.973183 


Bound (18) given in [22] 


0.967742 


0.993377 


NA 


0.360721 


0.0859599 


0.612245 



(iii) The sequence of Gegenbauer (or ultraspherical) polynomials {C^jJ^o is the one-parameter classical family 
obtained from the sequence of Jacobi polynomials {Pn'^}?^=o with a — P = X — ^. The sequence is 
orthogonal on [—1, 1] with respect to the weight function (1 — x^)^~^ for A > — ^ , n S N. The zeros of 
lie in the open interval (—1, 1) and are symmetric about the origin with a simple zero at the origin 
when n is odd. From [8, Theorem 2 (ii) d)] and Corollary 2.4, it follows that 



,,2 ^ 1 (2A + l)(2A + 3) 

^ n(n + 2A + 2) + (2A+l)(2A + 3)" ^ ' 



From (15), letting a = /3 = A — i, we obtain 



(2A + l)(2A + 5) 

~ 4n(n + 2A + 2) + (2A + 1)(2A + 5) ^ ^ 



Remark. The bounds obtained in (19) and (20) are particularly sharp when A is close to — ^. For n 
large, (19) and (20) yield wf > 1 - 0{^) and wi> 1- 0{^) respectively. 

In table 4 we compare (19) and (20) to the bound given by Szego (cf. [22, eqn. (6.2.13)]) 

^2>i_^A + l 
^- n + 2A + l' ^ ' 

noting that this lower bound also follows immediately from [8, Theorem 2 (ii) c)] by virtue of Corollary 
2.4. 



Table 4: Comparison of lower bounds for the largest zero of C^_^_i{x) for different values of n and A. 



Value of n and A 


n = 30, A = -0.49 


n = 3, A = -0.4956 


n = 50, A = 30.9 


n = 5, A = 60.49 


Largest zero of C^+i 


0.999978 


0.999267 


0.897416 


0.289296 


New bound (20) 


0.999978 


0.999267 


0.844369 


0.142782 


New bound (19) 


0.999978 


0.999266 


0.763561 


0.201482 


Bound (21) given in [22] 


0.999667 


0.998537 


0.66578 


0.198435 



4. Concluding remarks 

1. The bounds derived in this paper improve previously known upper (lower) bounds for the smallest (largest) 
zeros of Laguerrc, Jacobi and Gegenbauer polynomials, significantly so in the cases of both Jacobi and 
Gegcnbauer. Our results complement those obtained by other authors who derive upper (lower) bounds 
for the largest (smallest) zeros, notably, van Doorn in [23]; Dimitrov and Nikolov in [7] for the zeros of 
Laguerre, Jacobi and Gegenbauer polynomials; Ismail and Li in [15]; and Krasikov who gives uniform 
bounds for the extreme zeros of Laguerre polynomials in [16]. 

2. Various versions of (2) can be used to obtain sharper upper (lower) bounds for the smallest (largest) zeros 
of Laguerre, Jacobi and Gegenbauer polynomials by letting k = 3,4,... The calculations become more 
complicated as the degree of the coefficient polynomial Gk increases. 
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